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THE  MIXED  MODE  CRACK  PROBLEM  IN  A 
NONHOMOGENEOI  S  ELASTIC  MEDIUM^ 


by 


Noboru  Honda  and  F.  Erdogan 
Lehigh  L’niversity.  Bethlehem,  PA  18015 


ABSTRACT 

In  thij  paper  a  : onlioniogencc-ui  elastic  medium  containing  a  crack  arbitrarily  oriented 
with  respect  to  tlie  direction  of  property  gradient  is  considered.  The  problem  is  solved  under 
plane  strain  or  generalized  plane  stress  conditions.  This  is  a  highly  simplified  version  of  a  class 
of  physical  problems  that  may  arise  in  fracture  mechanics  studies  of  ceramic  coatings, 
metal/ceramic  composites  and  interfacial  zones  with  continuously  varying  volume  fractions  or 
graded  properties.  The  main  results  of  the  paper  are  the  calculated  modes  1  and  II  stress 
intensity  fetors.  Among  the  questions  studied  are  the  effects  of  the  material  nonhomogeneity 
constant,  the  crack  orientation,  the  loading  conditions,  and  the  Poisson’s  ratio  on  the  stress 
intensity  factors.  Briefly  discussed  are  also  the  stress  state  near  the  crack  tip  and  the  crack 
opening  displacement. 

INTRODUCTION 

Heterogeneous  materials  or  various  forms  of  composites  have  always  been  widely  used  in 
technological  applications.  These  materials  are  generally  designed  in  such  a  way  that  certain 
physical  bulk  properties  of  the  medium  are  optimized.  From  a  viewpoint  of  mechanical 
functioning,  in  practice  most  failure  processes  such  as,  for  example,  corrosion,  wear,  and 
fatigue  appear  to  be  surface  related.  Thus,  to  a  certain  extent  many  of  these  failures  can  be 
controlled  by  controlling  the  material  properties  near  and  at  the  surfaces.  In  some  cases  a 
relatively  simple  surface  treatment  would  be  sufficient  as,  for  example,  in  surface  hardening  to 
prevent  wear  or  in  introducing  residual  compressive  stresses  to  the  surface  to  prevent  fatigue 

^This  study  was  supported  by  .NSF  under  the  grant  MSM-8613611,  by  the  Office  of 
Naval  Research  under  the  contract  N00014-89-J-3188,  and  by  Sumitomo  Metal  Industries, 
Ltd. 

^Permanent  =».ddr?ss:  Sumitomo  Metal  Industries,  Ltd.,  Amagasaki,  Japan. 
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crack  initiation.  In  many  other  cases,  however,  a  protective  coating  with  a  more  resistant 
material  may  be  needed.  .A.s  some  examples  for  such  applications  one  may  mention  the 
thermal  barrier  ceramic  coating  of  combustion  chambers,  engine  blades  and  other  components, 
coating  of  machine  tools  against  wear,  coatings  used  for  protection  against  corrosion,  and  ion 
plating  of  parts  by  very  low  yield  materials  (such  as  gold  or  silver)  to  reduce  friction. 

In  designing  various  coatings,  in  addition  to  their  expected  physical  performance,  their 
mechanical  reliability  as  influenced  by  such  failure  related  factors  as  cracking  and  debonding 
must  also  be  taken  into  consideration.  To  a  large  extent  the  resistance  of  the  material  to  such 
failures  can  be  influenced  through  processing  techniques  (for  example,  by  controlling  the 
substrate  temperature  to  influence  the  residual  stresses,  and  by  controlling  the  degree  of 
mixing  in  plasma  spray  coating  or  in  ion  plating  to  influence  the  ductility  of  interfacial  zones) 

I  see  Batakis  anc  \  ogan.  19S5.  and  Houck.  1987).  It  is  at  this  point  that  influencing  material 
properties  through  tailoring  the  composition  as  well  as  controlling  the  processing  techniques 
may  be  a  practical  option.  It  has  been  shown  that  certain  strength  related  properties  of 
ceramic  coatings  can  be  improved  by  layering  the  interfacial  zone  going  from  metal  rich  to 
ceramic  rich  compositions.  For  example,  in  joining  tungsten  to  zirconia  by  introducing  four 
intermediate  layers  that  contain  80/20.  60/40.  40/60  and  20/80  percent  W/ZrOj,  respectively, 
it  wa'  shown  that  the  peak  value  of  the  residual  stress  becomes  approximately  one-sixth  of 
that  obtained  from  direct  \V-ZrO2  bonding  (Hirano,  et  al..  1988). 

The  next  logical  step  is.  of  course,  the  processing  of  fully  tailored  materials  and 
interfacial  zones  with  predetermined  continuously  varying  volume  fractions.  Advances  in 
powder  technology  and  surface  chemistry  in  recent  years  have  indeed  made  it  possible  to 
de\elop  such  nonhomogeneous  materials  (or  functionally  gradient  materials)  having  metal  and 
ceramic  constituents  (Hirano,  et  al.,  1988;  Hirano  and  Yamada,  1988).  A  different  type  of 
problem  in  which  one  would  have  to  consider  the  material  as  being  nonhomogeneous  would  be 
certain  thermal  stress  problems.  If  the  thermo-elastic  constants  are  significantly  dependent  on 
the  temperature  and  if  the  temperature  variation  in  the  medium  is  sufficiently  high,  then  for 
realistic  modelling  and  analysis  the  material  has  to  be  considered  as  being  nonhomogeneous. 

Fracture  toughness  and  fatigue  crack  growth  characterizations  of  the  nonhomogeneous 
materials  req’drp  the  solution  of  certain  standard  crack  problems.  The  mode  I  plane  strain 
problem  for  an  infinite  nonhomogeneous  medium  (that  is,  the  case  of  5=0  in  Fig.  1)  was 
considered  by  Delale  and  Erdogan  (1983).  who  showed  that  the  relative  dependence  of  the 
stress  intensity  factors  on  the  nonhomogeneity  constant  may  be  quite  significant.  Other  crack 
problems  involving  various  forms  of  ma’eria'  nonhomogeneity  were  considered  by  Kassir 
(  19721.  Dhaliwal  and  Singh  (1978),  GerasauUs  and  Srivastav  (1980),  Erdogan  (1985).  and 
Delale  and  Erdogan  (1988a.  1988b).  In  this  paper,  we  consider  the  generral  mixed  mode  plane 
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strain  problem  for  an  arbitrarily  oriented  crack  in  a  nonhomogeneous  medium.  In  previous 
studies  it  was  shown  that  the  effect  of  Poisson's  ratio,  v,  on  the  stres  intensity  factors  is  no; 
very  significant.  Thus,  in  this  study,  too.  u  is  assumed  to  be  constant. 

FOR.\IL  L.4TION  OF  THE  CR.ACK  PROBLEM 

Consider  the  plane  elasticity  problem  shown  in  Fig.  1  where  the  medium  contains  a 
finite  crack  on  y=0  plane  and  has  a  shear  modulus  fi  defined  by 
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fdo  and  6  being  material  constants.  By  observing  that  k=3-i/  for  plane  strain  and  k  =  (3- 
i/)/{l-ri/)  for  plane  stress,  and  hence 


A(x,yl  =  2^ 
^•1 
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(3) 


The  Xavier's  ecnations  for  the  elastic  medium  may  be  expressed  as 
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O'x-' 


dy 


^  =0  ’  (4a.b) 

where  u  and  v  are  respectively  the  x  and  y  ccomponents  of  the  displacement  vector.  After 
separating  the  solution  of  the  uncracked  medium  subjected  to  the  prescribed  external  loads, 
the  perturbation  problem  would  have  to  be  solved  under  the  following  boundary  conditions 
and  self-equilibrating  crack  surface  tractions  Pj  and  P2: 
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<Tyy(X.  +  0)  =  <Tyy(x.-0)  .  (Txy(X.+0)  =  <7xy(X,-0) 


(5a. b) 


v(x.+0)-v(x.-0)  =  0  .  |x|>a:  (Tyy(x.+0)  =  PjCx)  ,  |x|<a 


(6a.b ) 


u(x.-r0j-u(x.-0)  =  0.  |x|>a:  (Txy(x.-!-0)  =  P2(x)  ,  |x!<a 


(Ta.b; 


By  expressing  the  solution  of  (4)  as 


f  >: 

u(x.y)  =  ^  j  l'(y.a)e'“^^dQ 


,-fx.y)  =  ryz 


\'(y.Q  )e'‘^^dQ 


(8a, b) 


we  find 


■i  ‘i.y  4  n|3 

.  =  rnF(Q)e''  .  \(;y.Q)=.E  F:(Q)e  ■' 

j=l  J  J  j=l  J 


(9a.b) 


where  are  unknown  functions,  n^,....n4  are  the  roots  of 


n^  -f  2-'n^  +  i-2Q(Q-!-iJ)  +  "*•  Q'){-'^cx-\3 


+  Q2(Q2+2iQj-,j2  +  ^2  ^  0  ^ 


and  nij . are  given  by 


[i2o-)J(3-K)]nj+i(><-l)Q7 
•’  (x-l)nj^+(/c-l)7nj-(K+l)Q(a+i/^) 


,  j  = 


It  can  be  shown  that  the  characteristic  equation  (10)  may  be  expressed  as  follows: 


n^+7n-o(Q+iy?)]^  +  ^  (o7-i;?n)^  =  0  . 


From  (12)  it  then  follows  that 
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Since  u  and  v  must  vanish  for  x^+y"  —  x.  from  (9)  and  (13)  it  follows  that 
f'3!  o  )  =  I  4 !  o  '  =  0  .  y  >  0  . 

F  j  ( (1  =  F  2  (  u  '  =  0  .  y  <  0  . 

By  usina  tiio  Hu^)k^ '?  Law.  from  (S).  (9).  (15)  and  (16)  we  obtain 
,  r>^  f*l 

(Txxi  x.v  '  =  I  I  [-ium  (2^  +  A)-‘-njA]Fjexp(njy)exp(-iQx)dQ  . 

-  .  .  A-  j=f 
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[-ionijA  +  nj(2^i-r A)]Fjexp(njy)exp(-iQx)dQ  , 


f  +  1 

I  [njmj-iQ]Fjexp(njy)exp('iax)dQ  . 

-  X  j  =  ^ 


(IS) 


(19) 


where  £=1  for  y>0  and  £=3  for  y<0.  If  we  now  substitute  from  (18)  and  (19)  into  the 
homogeneous  conditions  (5).  we  find 

F3  =  RiFi  +  R2F2  ,  F4  =  R3F1  +  R4F2  (20a.b) 

where  the  known  functions  Ri(q) . R4(a)  are  given  by 

Rj(a)  =  {(  ni4-mj^)[(l-rx')t>i"4"h(3''f)Q^]+io(n4-n^)[l  +  «-(3-K)mj^m4]}/Ro  1 

Rjio)  =  {(m4-m2)[( l+'<)n2n4+(3-«)Q^]+io(n4-n2)[l+'c-(3-x)m2m4]}/RQ  . 
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R3(q)  =  -  {(mj-nij)'!  l-i-K)nin3— i3-K')o^j-i-iQ(n3-nj)[l-rK-i3-«)m]^m3]}/RQ  . 


RqIo)  =  ( 1114-1113  :’:  1 —  «  )n3n4-^(3-^•)Q^l+ia^n4-n3)[l-^^■-!3-^■j^l3m4]  .  ;21a-c 

The  two  remaining  unknown  functions  and  Fj  are  determined  from  the  mixed  boundary 
ccondirions  (6i  and  (7  1. 

THE  INTEGRAL  EQl  ATIONS 

To  reduce  the  probleii.  To  a  system  of  integral  equations  we  introduce  the  following  new 
unkiioun  functions: 


Jl 
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('22a.b) 

B>  substituting  for  u  and  \  from  the  results  found  in  the  previous  section  and  inverting  the 
Fourier  integrals,  from  '  it  niay  be  shown  that 

g^ft  'e"^'dt  =  -io  l-M'if^j-iof^jlFj  T  (fi2**af42)f2]/(‘^^)  • 

^  -  X 


gji  t  le'^'dt 


-i  j  { /I  I X )f2^ -f-iol 3-x )f33]F3  -f-  [*(l  +  x)f22 


^  io  ' 3-k  if32,r2  G  (  4Rq  1  , 


where 

fij  =  n3r"4mj(n4-nj '*n4m3mj(nj-n3)+njm3m4(n3-n4)  , 
^2j  ~  rii4mj(n4-nj '^m3mj(nj-n3)-fm3m4(n3-n4)  , 

^3j  ~  n3m3(  n4-nji^n4m4  nj-n3)-t-njmj(n3-n4)  . 
f4j  =  m3(n4-nj  l-m4(nj-n3)-rmj{n3-n4)  .  j=1.2. 


(23a.b) 


(24a-d) 


Equations  {23  >  and  f24)  give  the  unknown  functions  and  F2  in  terms  of  and  gj.  On  the 
other  hand,  b\  substituting  from  (IS)  and  (19)  into  (6b)  and  (7b)  we  find 

1  2  .  n-y-iox 

lim  ^  i-iom  A-rii|(2p't-A  )]F  (Q)e  do  =  Pilx;  .  |xi<a  . 

y  — +0  -■  j.x  j=l  '  ■'  ■'  •' 
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y-+0 


_L 

2t 


n.y-ioy 

Mi  njiiij-io  lFj(  o  )e  do  =  P2(x) 


ixi<a 


( 25a,b 


.^o\^  .  soK  ing  1 23 1  for  and  F2.  substituting  into  (25).  and  from  (22).  (6a)  and  (Fa; 
obser\  ing  that  g^i  t  '  =  0  for  t'>a.  we  obtain  the  following  integral  equations: 


lim  -  Im,!  x,\'.t  )g  !  t  )dt 

V--0  ...a  J  =  1  ^ 


t(k  — 1 ) 
-^o 


-3x 
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Pk(x)  .  |xi<a 


k=1.2  . 


(26) 


\v  hrr(> 


hj^j(x.y.t  '=  j"  Kj^jly.o.'r''^"^  ^^do  .  k=l,2;  j=l,2  ,  (27) 

and  the  known  functions  are  given  in  the  Appendix.  Since  the  conditions  gj(t)=0  rather 
than  (6a)  and  (6b!  arc  used  for  ;t|>a  in  deriving  (26),  the  integral  equations  must  be  solved 
under  the  following  single-valuedness  conditions: 

/■a 

g^ltldt  =  0  .  j  =  1.2  .  (28) 

-a 

In  order  to  determine  the  possible  singular  behavior  of  (26),  the  behavior  of  the  kernels 
h(jj.  (kj=1.2)  at  x  =  t  and  y=0  needs  to  be  examined.  For  this,  it  is  sufficient  to  determine  and 
separate  those  leading  terms  in  the  asymptotic  expansion  of  as  |o|— >oo  that  would  lead  to 
unbounded  integrals.  From  the  expressions  of  K^j  given  in  the  Appendix  it  can  be  shown  that 
in  the  asymptotic  expansions  for  o  — x  the  only  terms  that  would  give  unbounded  integrals 
are 


K.^^(y.a) 


1  M  -i®ly 

2i  a  ^ 


(29) 


Substituted  into  (27).  these  terms  give  Cauchy  type  kernels.  The  next  lower  order  terms  are 
of  the  form  l/loi  and  give  at  most  logarithmic  kernels,  log|t-x|  which  are  square  integrable  and 
may  be  treated  as  Fredholm  kernels.  By  adding  and  subtracting  the  asymptotic  values  given 
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bs’  (29)  to  and  from  in  (27).  and  by  evaluating  the  integrals  involving  the  leading  terms. 
(26)  may  be  modified  as  follows: 


I  kiiix.t)gi(t)+ki2(x.t)g2(t)]dt  =  Pi(x)  .  |x|<a  . 

^  k2i(x.t)gi(t)-k22(x.t)g2(t)]dt  =  P2(>«)  •  Ix|<a  . 

where 

k^jlx.t  i  =  hj^j(x.O.t)  .  k22fx.t)  =  h22(x.0.t)  . 
k|2(x.t  )  =  I  ■  K ^ (O.o  )’e  ^  • 

J  -  X 

k2i(x.t  .=  [K2i(  O.a  )-K2,_^^(0.a  )je'‘^^’^'^^dQ  . 


(30a.b) 


(31a-d) 


SOLUTION  .\ND  THE  STRESS  INTE.NSITY  FACTORS 

To  solve  the  system  of  integral  equations  (30)  we  first  define  the  following  normalized 
quantities: 

s  =  t/a.  r=x/a.  g,( t )  =  f,(s').  k,j(x,t)=:Ly(r,s).  Pi(x)=Pj(ar),  (ij  =  1.2).  (32) 


Then,  for  example.  (30a)  may  be  expressed  as 


1  f'  ^  2 


jli  L,j(r.s)fj(s)]ds  =  pi(ar),  |r|<l 


(33a) 


.2.-1/2 


Noting  that  the  fundamental  solution  of  (33a)  is  (1-s^)  ,  the  unknown  functions  fj^  and  f2 

may  now  be  expressed  as 


fj(s)  =  I  .  f2(s)  =  I  , 


n=0  n-s^ 


(34a, b) 


where  Tn  is  the  Chebyshev  polynomial  of  the  first  kind.  From  (28),  (34)  and  the 
orthogonality  conditions  of  Tn(s)  it  may  be  seen  that 

Ao  =  0  .  Bo  =  0  .  (35) 
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kj(-a)  (Tyy(x.O)  .  k2(-a)  ^  •j2(-x-a)  <Txy(x,0)  .  (38a-d) 

To  evaluate  kj  and  k2  "e  observe  that  equations  (25),  or  (30)  provide  the  expressions  for 
(7yy(x.0)  and  (r^yl^-O)  that  are  valid  for  |x|>a  as  well  as  |x|<a.  Thus,  by  using  (36),  and 
replacing  Pi(x)  by  cr yy(x.O)  (for  |x|>a).  from  (33)  it  may  easily  be  shown  that 
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HESL  L  I> 

1  lie  main  results  of  thi-  study  are  tlie  stress  intensity  factors  that  are  obtained  from 
(40  ;  and  i41)  after  solving  equations  (37),  and  are  given  in  Tables  1-7.  Referring  to  Fig.  1  we 
observe  that  the  half  crari.  length  a  is  the  only  length  parameter  in  the  problem  which  may  be 
eliminated  by  suitably  normalizing  all  relevant  quantities  (see,  for  example,  Eq.  32).  The 
exponent  i  in  pi,Xj)  =  poexp(i5x^)  is  the  measure  of  material  nonhomogeneity  and  appears  as 
the  nondimensional  constant  a.6  in  the  results  given.  The  angle  9  defines  the  orientation  of  the 
crack  with  respect  to  the  direction  x^  of  the  material  property  gradient.  As  physically 
expected,  for  6=0  the  integral  equations  uncouple  and  the  crack  problems  corresponding  to 
mr  ues  I  and  II  loading  conditions  can  be  solved  separately.  For  9^0  the  integral  equations  are 
always  coupled  regardless  of  the  loading  conditions.  Tables  show  the  stress  intensity  factors 
k,  and  k2  in  normalized  form.  The  basic  definition  of  and  k2  are  given  bv  (38).  The 
normalized  stress  intensity  factors  are  defined  by 

k,<  =.a)  =  k,(4:a)/(PN^)  (42) 
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where  P  represents,  in  each  case,  the  amplitude  of  the  crack  surface  traction.  Referring  to 
(bb'  and  (7b;.  the  following  tractions  were  used  in  the  examples  given  by  the  Tables: 

P2(x)  =  0  .  Pitx)  =  (-(Tq.  -c^jx/a  .  -{TjX^/a^  .  -^3X^/3^;  .  (43; 

Pji X )  =  0  .  P2' X;  =  (-Tq  .  -  r^x/a  .  -rjX^/a^  .  -r3X^/a^J  .  (44  j 

For  crack  surface  tractions  that  are  reasonably  smooth  functions  of  x.  the  stress  intensity 
factors  may  be  approximated  by  a  suitable  superposition  of  the  results  given  by  (43)  and  (44  j. 
It  should  be  pointed  out  that  in  nonhomogeneous  materials  the  external  loading  conditions  in 
mar.>  cases  ma>'  be  modelled  by  prescribing  displacements  rather  than  tractions.  The  two 
practical  cases  that  may  be  considered  are  the  “uniform  strain"  and  "bending"  applied  to  the 
mediu.m  a\sa>  from  the  crack  region.  It  will  be  assumed  that  the  displacements  are  prescribed 
in  the  planes  parallel  to  the  direction  of  the  property  gradient,  x^  in  such  a  way  that 


-12- 


)  =  corresponding  crack  surface  tractions  may  be  shown  to  be 


'yyl 


x.O)  =  -(Tosin^f/  .  cxy'x.O)  =  -iTosin^cos^ 


(-19; 


whicli.  for  a  giveii  crack  orientation  6.  are  constant.  Thus,  the  solution  for  the  loading 
conditions  (49)  may  be  obtained  by  a  proper  superposition  of  the  results  found  by  using  (43; 
and  ( 44  i. 

Table  1  shows  the  effect  of  aif  on  the  stress  intensity  factors  for  the  two  limiting  crack 
orientations.  ^  =  0  and  ^  =  t/2  under  a  uniform  crack  surface  pressure  cto-  For  ^=0  the  crack 
results  obtained  are  es'eutially  the  same  as  that  found  by  Delale  and  Erdogan  (1983;.  For 
0=-/2.  x=0  is  a  plane  of  s>nimotry.  the  crack  problem  is  one  of  mixed  mode,  and  the  mode  II 
component  becomes  quite  significant  for  larger  values  of  aS.  For  the  uniformly  pressurized 
crack  more  detailed  results  are  shown  in  Table  2  where  the  angle  6  is  varied  between  0  and 
Tr/2.  The  results  need  to  be  examined  rather  carefully  if  one  is  interested  in  the  fracture 
initiation  at  the  crack  tip>.  For  example,  one  surprising  result  shown  by  Table  2  is  that  the 
maximum  values  of  k^  and  kj  do  not  generally  correspond  to  the  limiting  crack  orientations 
<1=0  and  d~rri'2.  This  may  be  seen  somewhat  more  clearly  from  Fig.  2  which  shows  the 
variation  of  the  stress  intensity  factors  with  the  angle  9  for  a6=l.  Also,  the  values  of  6 
corresponding  to  maximum  stress  intensity  factors  seem  to  depend  on  a6. 

For  a  quick  assessment  of  a  possible  crack  growth  initiation,  it  is  generally  sufficient  to 
examine  the  aniplitude  and  the  direction  of  the  maximum  cleavage  stress  at  the  crack  tips 
^see  Erdogan  and  Sili.  19Gj)-  Given  the  stress  intensity  factors  kj  and  k2,  the  "cleavage" 
stress  in  the  small  neighborhood  of  the  crack  tip  may  be  expressed  as 


(T  ^  (r.c) 
O 


"k.  CO'* 


o  3 

O  ■  o 


k2sind]cos 


0 


(50) 


where  (r,<p)  are  the  polar  coordinates  at  the  crack  tip  and  0  is  measured  from  the  x  axis. 
Generally  the  hypothesis  is  that  the  crack  initiation  would  be  radial  in  a  direction  0* 
perpendicular  to  the  local  maximum  cleavage  stress  obtained  from 


A 

ao 


[kjsino*  +  k2(3cos<>*-l)]  =  0  . 


cr  (  r.O*  !  >  0 
00 


(51a.b) 
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and  at  a  value  of  the  load  level  determined  bv 


(r.<5 

oo 


=  A 


o 

O 


3  k 
0 


0 

sino  )cos  = 


>l2Tr 


(521 


where  is  the  (local  i  critical  stress  intensity  factor  of  the  medium.  In  studying  the  fracture 
initiation  from  an  e.xisting  flaw  in  brittle  and  quasi-brittle  non  homogeneous  materials,  two 
additional  points  need  to  be  made.  The  first  concerns  the  angular  distribution  of  the 
asymptotic  stress  state  for  small  values  of  r.  It  appears  that  these  expressions  (such  as  (50)) 
for  a  nonhomogeneous  medium  are  identical  to  that  of  a  homogeneous  medium  provided  near 
and  at  the  crack  tip  the  elastic  properties  of  the  medium  are  continuous  (but  not  necessarily 
difleret.tiablej  functions  of  the  space  coordinates  (see  Delale  and  Erdogan,  1988a).  The  second 
point  i'  that,  e\en  though  the  material  is  assumed  to  be  isotropic  with  respect  to  its  fracture 
resistance,  the  resistance  parameter  is  expected  to  be  a  function  of  the  space  variables. 

may  be  determined  from  standard  fracture  toughness  experiments  by  using  a  series  of 
homogeneous  specimens  covering  the  complete  range  of  material  composition  for  the 
nonhomogeneous  medium  under  consideration. 

For  a  fixed  value  of  a^=l  Table  3  shows  the  effect  of  loading  conditions  and  the  crack 
orientation  on  the  normalized  stress  intensity  factors.  It  may  be  observed  that  under  certain 
individual  loading  conditions  the  mode  I  stress  intensity  factor  could  be  negative,  implying 
crack  closure.  Such  solutions  are,  of  course,  not  valid  and  can  only  be  useful  in  a  superposition 
that  gives  a  positive  resultant  k^. 

In  the  solution  given  effect  of  the  variation  of  the  Poisson's  ratio  v  is  assumed  to  be 
negligible.  .Analytically,  it  is  difficult  to  verify  this  assumption.  However,  one  can  solve  the 
problem  for  various  different  values  of  v  and  compare  the  results.  This  is  shown  in  Tables  4 
and  5.  Table  4  shows  the  effect  of  v  on  the  stress  intensity  factors  for  aj=0.25,  tf=x/2,  and 
for  uniform  crack  surface  tractions  -<Tq  and  -Tq.  The  maximum  difference  observed  was  2%  in 
kj  under  uniform  crack  surface  pressure  (Tq.  The  difference  for  other  loading  conditions  (which 
are  not  all  shown  in  the  table)  was  somewhat  smaller.  For  9—0  similar  results  were  found  by 
Delale  and  Erdogan.  1983.  The  effect  of  v  is,  however,  more  significant  for  greater  values  of 
a<$.  For  example,  the  same  calculations  as  Table  4  were  repeated  for  aj=2.5  and  the  results 
are  shown  in  Table  5.  It  may  be  seen  that  the  difference  in  kj  obtained  from  0.05<i/<0.45 
may  be  as  high  as  15%.  On  the  other  hand,  within  a  more  practical  range  of  i/,  namely  for 
0.2<i/<0.35.  the  difference  is  less  than  6%. 

Some  additional  results  for  displacement  loading  given  by  (45)  are  shown  in  Tables  6 
and  7.  The  tables  show  the  effect  of  the  non  homogeneity  parameter  ai  and  the  crack 
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orientation  9  on  the  normalized  stress  intensity  factors. 

.\fter  solving  the  integral  equations  and  observing  that 


r  = .  1 V 

1  J  i.s2  n  ^ 


the  relative  crack  surface  opening  may  easily  be  obtained  from  (22),  (32),  and  (34)  as  follows: 


(■X  I -  3^. 

u(  x.+0)-u(x.-0)  =  gift  )dt  =  -  E  ^  AnL'j^.2(x/a)  , 

J  -a  1 

- -  X 

V,  x.-r0,]-v(x.-0  :■  =  g2it)dt  =  -xa^-x^  E  i  BnU  ,(x/a)  . 

J  -a  1 


(54a.b) 


Figures  3  and  4  show  the  crack  opening  in  y  direction  for  6=0.5  and  2.5,  respectively,  where 
the  loading  is  crack  surface  pressure  ctq  and  the  normalized  displacements  given  in  the  figures 
are  defined  by 

^  [''IX-Ol-vlx.-O;:  .  Vo  =  (To  ^  .  (55) 

The  three  displacements  shown  in  each  figure  correspond  to  S=0  and  9=^ir/2  with  696O.  and  to 
a  homogeneous  material  (6=0).  It  may  be  seen  that  for  large  values  of  6  the  crack  opening 
displacements  in  nonhomogeneous  materials  can  be  significantly  greater  than  the  corresponding 
homogeneous  values. 
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APPENDIX 


Expressions  of  the  functions  K,j(y.Q).  (ij  =  l,2} 


i\, , !  y .a  I  =  e  '  ■  { 


*\  ,  i  3-K)Qni,  — i(  l-t-x-in, 

—  '  J  ^  * 


(3-x  lQm-,-‘-i(  l-rKin-,  ,  ,  ,  n^v, 

i )  , 


f  I jam, ^  ^  ^ 

=p  '  i  - ; — i— - *  f 

“  4a  ‘  i  /-.Q  ^ 


- — 2  (l^x-)(f,,-.af4i)e  ■}. 


, .  ,  -v  %•,  Q  -^111 ,  m .  ,  ,  ,  n,  V 

Kji  .v.ci  )=  e '•  {  — i  ;;  1-K!f32-1Q(3-K)f22]e 


a-i-in2ni 

■4o  -o 


-  .•■  1-X  Ifji 


ia(3-K)f2i]e”^''  }  . 


1\2'.'  >  .u  .1  =  e  ■  { 


V  ,  0-111,  m,  ,  .  r  ■  "iV 


Q  ^iiioni 


4q-* 


2 — 2  ii-,.«)(f  .iQf^^)e  } 


-c  =  '  mi-ni2){ni3-m4')(njn2-^n3n4)  +  (m3-m4)(m2-m3)(n2n3+njn4) 


-  irni-m3)im2-m4)(njn3T-n2n4) 


where  y.  m|.  Uj,  and  fj,,  are  given  by  (2),  (11).  (13),  and  (24).  respectively. 
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Table  1.  The  effect  of  the  nonhomogeneity  constant  a:  on  the  stress 
intensity  factors;  '.=0.3,  p  (x)=-r  ,  p  (x)=0,  k.(+a)  = 
k.(+a)/r  va  .  1  o  2  i 

1  o 


-  = 

-/2  (k^(-a)  = 

=  k^(a). 

k^C-a)  =  -k^ 

(a)) 

a  • 

0.1 

0.25 

0.5 

1.0 

2.5 

5.0 

(a) 

1 

1.008 

1.036 

1.101 

1.258 

1.808 

2.869 

k^  (a) 

.  0.026 

0.065 

0.129 

0.263 

0.697 

1.567 

-  =  0  *  0’ 

o 

H 

1 

k^(a'i 

1.023 

1.055 

1,103 

1.189 

1.382 

k^(-a) 

■  0.975 

1 

0,936 

0.871 

0.757 

0.536 
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intensity  factors;  =  0.  J,  |i  (x)  =  -o  ,  |>,,(x)  ==  0,  k.(f.i)  =  k.i+.i)/' 


19 


065  -0.099  -0.138  -0.181  -0.223 


Table  1.  The  effect  of  0  and  the  lo.uiin>’  conditions  on  the  stress  int 
a.S  =  1.0,  k  (+a)  =  k  (+a)/ (c  .  ,  I  .  )  . 
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k»(-a)  0  -0.072  -0.143  -0.208  -0.264  -0.305  -0.321  -0.303  -0.244  -0.141 
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Table  1.  The  effect  of  .  on  the  stress  intensity  factors;  a‘>0.25,  -=-/2, 

=  k.(+a)/(r  )ta  . 

1  1  0  0 


0.05  0.1 

0.15 

0.2 

0.25 

0.3 

0.35 

0.4 

0.45 

p^vx)  =  -- 

P-,(X)  =  0,  (k 

= 

k^(a),  k 

,(-3)  = 

-k,  (a) ) 

k^(a) 

;  1.026  1.025 

1.030 

1.032 

1.03- 

1.036 

1.039 

1.042 

1.046 

rC  T  (a) 

;  0.065  0.065 

0.065 

0.065 

0.065 

0.065 

0.065 

0.065 

0.065 

F-,  (x'l  =  0, 

p  ^  (  X  1 

=  (k 

c 

./-a)  = 

-k^(a) , 

k.iC-a') 

=  k,(a) ) 

k^  (a  ■) 

i-0.063  -0.063 

-0.063 

-0.063 

-0.062 

-0.062 

-0.062 

-0.062 

-0.062 

( a ) 

0.995  0.995 

0.996 

0.996 

0.997 

0.998 

0.999 

1.000 

1.001 

Table 

5.  The  effect  of 

and 

the  loadint  conditions  on 

the  stress  intensity 

factors;  a  = 

-  T.  X .  ("sa  1  -  k  . 

(+a)/(r. 

,  t  . )  ‘  a 

1 

1 

3 

3 

0.05  0.1 

0.15 

0.2 

0.25 

0.3 

0.35 

0.4 

0.45 

p,  (xl  =  -' 

c’  ^2^-- 

)  =  0,  (k 

j^(-a)  = 

kj^(a),  k 

2(-a)  = 

k2(a) ) 

k^(a) 

'  1.67-  1.696 

1.720 

1,747 

1.776 

1.808 

1.843 

1.882 

1.927 

k-,(a) 

0.6S5  0.6S" 

G.6S9 

0.691 

0.694 

0.697 

0.700 

0.703 

0.705 

r  ,  V  >; ,)  =  -  - 

,x  a,  Pt(x')  =  0, 

( (—a) 

=-k^(a) , 

k^C-a) 

=k,(a) 

kj_(a') 

0.633  0.63S 

0.6-3 

0.649 

0.656 

0.663 

0.671 

0.680 

0.690 

ko  (a) 

0.082  0.083 

0.084 

0.085 

0.086 

0.087 

0.088 

0.090 

0.091 

Pl(x)  =  -: 

,x--/a- , 

0,  (k^{- 

a)=kj^(a) 

,  k2(-a)=-k2(a) ) 

tc^  (a ) 

0. 66S  0.674 

1 

0.680 

0.687 

0.694 

0.703 

0.712 

0.722 

0.733 

(a) 

:  O.l'l  0.174 

0.175 

0.175 

0.176 

0.176 

0.177 

0.177 

0.178 

1  Pj(x)  =  0, 

=  (k 

^(-a)*-k 

:^(a),  k2 

(-a)*k2 

(-a)) 

k^(a) 

1 

i-0.405  -0.404 

-0.402 

-0.400 

-0.399 

-0.397 

-0.395 

-0.392 

-0.390 

k,(a) 

!  0.989  0.996 

1 

1.003 

1.011 

1.019 

1.028 

1.038 

1.048 

1.060 

^  Pj^(x)  =  0, 

=  -Tj^x/a, 

(kj^(-a) 

*kj^(a) , 

it^f-a)* 

-k2^a>) 

kj^(a) 

-0.020  -0.020 

1 

-0.020 

-0.020 

-0.021 

-0.021 

-0.021 

-0.022 

-0.022 

k-,(a) 

.  0.527  0.528 

0.529 

0.531 

0.532 

0.533 

0.535 

0.537 

0.539 

Pj^(x)  =  0, 

Pt(x) 

=  --^x'/a 

%  (k^(- 

■a)=-kj^(a),  ic^C- 

a)*k2(a) ) 

K^(a) 

-0.106  -0.105 

-0.105 

-0.104 

-0.104 

-0.103 

-0.103 

-0.102 

-0.102 

k-,(3) 

1  0.499  0.501 

0.503 

0.505 

0.508 

0.510 

0.513 

0.516 

0.519 
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Table  6. 

Normalized  stress 

intensity  factors 

for  "uniform 

strain" 

'  V  1  V 

=  away  from  the  crack 

region;  ■  =0. ; 

3, 

K  '  = 
0 

E  ►  a  . 

0  0 

c  a 

/  - 

k^U)/K 

1  0 

k,  (-a)/K 

1  o 

k.(a)/K 

2  0 

k,(-a; /K 

1  0 

1 

0 

'  1.196 

0.825 

0 

0 

i 

1 

0.1 

1.081 

0.750 

-0.321 

-0.254 

0.25 

0.2 

0.781 

0.548 

-0.514 

-0.422 

0.3 

;  0.414 

0.290 

-0.504 

-0.437 

0.4 

o.i;i 

0.075 

-0.304 

-0.282 

0.  5 

0 

0 

0 

0 

0 

1.424 

0.674 

0 

0 

' 

0.1 

1.285 

0.617 

-0.344 

-0.213 

0.2 

0.925 

0.460 

-0.548 

-0.365 

0.5 

t 

0.3 

0.490 

0.247 

-0.532 

-0.397 

0.4 

0.146 

0.059 

-0.314 

-0.269 

0.  5 

0 

0 

0 

0 

0 

6,317 

0.115 

0 

0 

0.1 

5.376 

0.117 

-0.867 

-0.037 

0.2 

3.315 

0.115 

-1.155 

-0.090 

2.5 

0.3 

1.441 

0.082 

-0.900 

-0.158 

0.- 

0.369 

0.004 

-0.429 

-0.179 

0.5 

i  0 

0 

0 

0 
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Table  7. 


Norroalized  stress  intensity  factors  for  "bending" 
the  crack  reaion;  ;=0.3,  c  (x,  ,+^)  =  t-iX,,  K,  = 

y,y,  1»  1  1 


away  from 
EtE  va  . 

1  o 


:  a 

1 

!  -/- 

k^(a)/K^ 

k^(-a)/Kj^ 

k  (a) /K, 

L  1. 

k2(-a)/K^ 

;  0 

0.637 

-0.391 

0 

0 

0.1 

0.5-2 

-0.340 

-0.174 

0.111 

0.2 

0.  il- 

-0.217 

-0.232 

0.158 

0.25 

0.3 

O.llS 

-0.088 

-0.160 

0.121 

0.4 

o.oie 

-0.014 

-0.049 

0.042 

'  0.5 

0 

0 

0 

0 

0 

o.so? 

-0.304 

0 

0 

o.i 

0.6S3 

-0.268 

-0.214 

0.087 

O.i 

0.2 

0.397 

-0.178 

-0.278 

0.128 

0.3 

1 

0.139 

-0.076 

-0.184 

0.105 

■  0.4 

0.018 

-0.013 

-0.053 

0.039 

0.5 

0 

0 

0 

0 

0 

4.982 

1 

-0.039 

0 

0 

!  0.1 

3.929 

-0.039 

-0.898 

0.012 

!  0.2 

1.880 

-0.037 

-0.968 

0.024 

2.5 

'  0.3 

0.486 

-0.028 

-0.480 

0.033 

0.4 

G.0-3 

-0.009 

-0.094 

0.023 

0. 5 

0 

0 

0 

0 
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1 


in  the  nonhor.ogeneous  r.ec 


*) 
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